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ABSTRACT 

A  study  is  made  of  the  hydrology  of  watersheds  less  than  200  square 
miles  in  area  located  in  the  State  of  Indiana  for  which  flows  are  recorded 
by  the  U.S.G.S.  A  statistical  frequency  analysis  of  the  peak  flows  was 
made  by  means  of  Gumbel's  extreme  value  theory.   A  geomorphological  study, 
involving  the  quantitative  determination  of  five  watershed  characteristics, 
combined  with  the  use  of  multiple  correlation  techniques  serves  to  establish 
a  formula  for  estimating  the  peak  discharge  of  ungaged  watersheds  in 
Indiana. 


SYNOPSIS 

The  purpose  of  this  research  is  to  study  the  hydrology  of  watersheds 
less  than  200  souare  miles  in  area  throughout  the  State  of  Indiana*and  to 
establish  methods  of  estimating  the  peak  discharge  from  these  watersheds. 
The  study  will  include  four  parts:  a  statistical  analysis  of  existing 
data,  a  geomorphological  study  of  selected  watersheds,  a  study  of  small 
watershed  storm  hydrographs,  and  an  analysis  of  overland  flow.  The 
present  paper  summarizes  the  work  done  on  the  first  two  parts. 

The  statistical  analysis  has  been  made  by  means  of  the  extreme 
value  method.  The  25-year  instantaneous  peak  flow  has  been  determined  for 
all  the  gaging  stations  with  a  drainage  area  of  less  than  200  square  miles 
for  which  the  U.S.G.S.  publishes  records.   For  those  watersheds  for  which 
sufficient  topographic  data  were  available,  five  geomorphological  parameters 
were  calculated:  the  watershed  area,  the  mean  relief,  the  main  stream  slope, 
the  stream  density,  and  the  watershed  shape  factor.  A  multiple  correlation 
was  made  between  these  variables  and  the  25-year  instantaneous  peak  flow. 
A  correlation  chart  is  given  to  obtain  graphically  the  peak  flow  for  small 
watersheds  from  these  five  geomorphological  parameters.  A  simple  formula 
containing  only  three  geomorphological  parameters  was  introduced  as  a 
first  approximation  of  flood  peak  discharge  determination.  An  expression 
is  given  to  calculate  the  peak  discharge  for  any  frequency  from  the  25-year 
peak  discharge.  Similar  analyses  could  be  done  for  other  regions. 


INTRODUCTION 

The  determination  of  the  required  waterway  area  of  a  bridge  or 
the  selection  of  the  size  of  a  culvert  are  problems  which  require  an 
accurate  estimate  of  the  peak  flood  discharge  that  will  pass  through 
the  structure.  The  determination  of  this  discharge  is  more  difficult 
for  small  watersheds  because  the  majority  of  them  are  ungaged.  There 
are  very  few  gaged  small  watersheds  on  which  to  base  an  estimate.  Par- 
ticularly in  the  state  of  Indiana,  there  is  very  little  information 
on  watersheds  less  than  200  square  miles.  There  are  only  twelve  water- 
sheds of  less  than  100  square  miles,  and  seventeen  watersheds  with  an 
area  between  100  and  200  square  miles  for  which  the  U.S.G.S.  is  currently 
reporting  flows,  a.  research  program  was  initiated  at  Purdue  University 
t»  obtain  reliable  methods  for  estimating  the  peak  discharge  for  a  safe 
and  economic  design  of  highway  drainage  structures,  serving  watersheds  in 
Indiana  of  less  than  200  square  miles  but  larger  than  20  square  miles. 

The  existing  methods  are  empirical  and  fail  to  take  into  account  the 

(3)  U£ 

factors  upon  which  the  runoff  depends.  Kinnison    in  1946  and  Chov* 

in  1962  have  given  a  complete  list  of  empirical  formulas  which  include 
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Kinnison,  H.  B. ,   "Flood-Flow  Formulas,"  Jour.  Boston  Soc.  Civil  Engr., 
V.  33,  pp.  1-19,  January,  1946. 

4 
Chow,  V.  T.,  "Hydrologic  Determination  of  waterway  Area  for  the  Design 
•f  Drainage  Structures  in  Small  Drainage  Basins,"  Engineering  Experi- 
ment Station  Bulletin  No.  462,  University  of  Illinois,  1962. 


the  watershed  characteristics.   The  most  frequently  used  formulas  are 

(5)  (6) 

those  of  Talbot   published  in  1887  and  keyer    in  1879  and  the 

(7) 
Kational  formula  originally  derived  by  kulvany    in  1857.  Talbot's 

formula  was  originally  intended  for  locations  in  Illinois.   It  estimates 

the  waterway  area  from  the  watershed  area.  The  formula  is: 

A  =  CM  ?/4  (1) 

where  A  is  the  required  waterway  area  in  square  feet,  h   is  the  watershed 

area  in  acres,  and  C  .is  a  coefficient  varying  between  1/5  and  1  depending 

on  the  slope  and  character  of  the  watershed.  The  selection  of  the  coefficient 

depends,  among  other  things,  on  the  experience  of  the  designer.  Due  to 

the  various  factors  that  affect  the  runoff  other  than  the  watershed  area,  the 

value  of  the  coefficient  C  cannot  be  accurately  determined  to  represent  all 

the  watershed  characteristics.   Talbot's  formula  is  unsatisfactory  for  a  safe 

design  of  a  hydraulic  structure. 

(a)  * 

Yule    in  1950  developed  a  similar  formula  for  locations  in  Indiana.  The 
required  waterway  area  is  expressed  as  a  function  of  the  2/3  power  of  the 
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Talbot,  a.  N. }    "The  Determination  of  .aterway  for  Bridges  and  Culverts," 

3e'lccteJ  P-:;-crs  of  the  Civil  jn  inooy'»»Club,  Tochno  m-h  ?Jo.2, 
University  of  Illinois,  1G87-C  .  ;  .  14-22, 
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1'ieyer   The  formula  was  first  published  in  a  paper  read  by  Cleemann  before 

the  Engineers '  Club  of  Philadelphia  in  1879* 
Cleemann,  T.  k. ,,  "Railroad  engineers'  Practice,  Discussion  of  Formulas," 

Proceedings,  Engineers  Club  of  Philadelphia,  V.  1,  pc.  146,  1879. 

7 
kulvany,  T.  J.,  "On  the  Use  of  Self  Kegistering  Rain  and  Flood  Gauges 

in  kaking  Observations  of  the  Relations  of  Rainfall  and  of  Flood 
Discharges  in  a  Given  Catchment,"  Transactions  of  the  Institution 
of  Civil  Engineers  of  Ireland,  V.  4,  Part  II,  pp.  18,  1857. 
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Yule,  B.  B.,  "Bridge  waterway-Area  Formula  Developed  for  Indiana," 
Civil  Eng.  V.  20,  pp.  26,  73,  Oct.  1950. 


drainage  area.  It  is 


A  =  Ca2/3  (2) 


••'here  A  is  the  waterway  area  in  snuare  feet,  a  is  drainage  area  in  square 
miles,  and  C  is  a  general  coefficient  which  varies  with  the  watershed 
topography:  C  is  0.3  to  0.7  for  flit  land,  0.7  to  1.3  for  rolling  land, 
and  1.3  to  2.2  for  hilly  land. 

Benson'-  °'  in  1959  found  that  the  mainstream  slope  is  next  in  importance 
to  drainage  area  among  the  factors  which  affect  the  runoff.  The  slope  con- 
sidered was  for  that  part  of  the  mainstream  located  between  85  to  10  percent 
of  the  total  distance  above  the  gaging  point.  The  following  empirical  formula 
was  found  for  the  New  England  region: 

Q.  =  a.Ab  sc  (3) 

where  Q  is  peak  discharge  in  cfs,  A  is  drainage  area  in  snuare  miles,  S 
is  th«  S5o-10^  slope  of  the  mainstream,  and  a,  b,  c,  are  the  regression 
coefficients  varying  with  the  recurrence  interval  of  flood. 

The  factors  affecting  the  runoff  may  be  grouped  in  three  categories: 
the  storm  characteristics,  the  geomorphological  and  the  geological  character- 
istics of  the  watersheds.  The  dependence  of  the  runoff  on  the  geomorphology 
of  the  watershed  is  analyzed  in  this  paper.  The  effect  of  the  storm  char- 
acteristics and  of  the  soil  types  are  considered  in  the  following  paper  en- 
titled "Design  Hydrographs  for  Small  '  atersheds  in  Indiana".  The  validy 
of  the  unit  hydrograph  theory  will  ultimately  be  analyzed  by  a  theoretical 
study  of  overland  flow  hvdrodynamics . 


(9)  Benson,  M.  A.,  "Channel-Slope  Factor  in  Flood-Freouency  Analysis," 
Jour.  Hydr.  Div.  A.S.C.E.,  April  1959. 


STATISTICAL  ANALYSIS 

The  following  two  methods  give  a  linear  relationship  between  the  ob- 
served variable,  i.e.  the  discharge,  and  the  recurrence  interval  or  a 
function  of  the  recurrence  interval. 

Geyer,(10'  in  1940,  derived  a  mathematical  expression  relating  the 
flood  magnitude  to  the  exceedance  interval 

Y  =  LR"tk  (4) 

where  Y  is  the  flood  magnitude  in  sec-ft.  that  has  an  exceedance  interval 
of  t  years,  and  L,  R,  Y   are  constants  for  a  particular  stream.  Taking 
logarithms  twice  on  both  sides  of  Eq,  4,  one  obtains 

log  (log  L  -  log  I)  -  K  log  t  +  log  log  R  (5) 

where  L  is  then  determined  by  trial  and  error  so  that  Eq.  5  represents  a 
straight  line  when  log  (log  L  -  log  Y)  is  plotted  against  log  t. 

Gumbel,  (^J-)  A 12),  (13)  ±n   1941 f  developed  the  extreme  value  theory 
which  is  used  to  analyze  the  observed  extremes  and  to  forecast  further 
extremes.  The  theory  states  that  the  probability  c'Xx)  for  the  discharge 


(10)  G-eyer,  J.  C,  "New  Curve-fitting  Method  for  Analysis  of  Flood  Records," 

Trans.  Amer.  Geophys.  Union,  1940,  II,  pp  660-666. 

(11)  Gumbel,  E.  J.,  "Statistics  of  Extremes,"  New  York,  Columbia  University 

Press,  1958. 

(12)  U.  S.  Dept.  of  Commerce,  "Probability  Tables  for  Analysis  of  Extreme- 

Value  Data,"  National  Bureau  of  Standards,  Applied  Mathematics 
Series  22,   1953. 

(13)  U.  S.  Dept.  of  Commerce,  "Statistical  Theory  of  Extreme  Values  and 

Some  Practical  Applications,"  National  Bureau  of  Standards, 
Applied  Mathematics  Series  33,  1954. 


x  to  be  the  largest  among  n  independent  observations  is  given  by 

(^(x)  =  e-e-y  =  exp  (-e~y) 

(6) 
as  n  =>  oO 

where  e  is  the  base  of  Napierian  logarithms  and  y,  termed  the  reduced 

largest  value,  is  given  by 

where  c*  ,  u  are  two  extreme  parameters ,  u  is  a  certain  expected  largest 

value  having  the  return  period  n,  and  its  probability  K^O  is  1-  =J 

n        n 

cy  is  defined  as  nf(un)  where  f(un)  is  the  initial  distribution  given  by 
f(un)  =  F'Cu^).  A  probability  paper  designed  for  extreme  value  was  pro- 
posed by  Powell. ^^J  The  observed  variate  x  is  traced  on  the  ordinate,  the 
largest  value  y   is  traced  on  the  abscissa,  both  in  linear  scales.  The  value 
of  the  probabilities  fTj(x)  is  piven  in  En.  6  or  by 


<I>(x)  = 


m 


n  +  1  (8) 

Eo,  8  gives  the  "plotting  position"  v/here  m  is  the  rank  of  the  yearly  max- 
imum in.  Increasing  order,  and  n  is  the  number  of  years  of  observation. 
The  return  period  is  given  by 

T=  1  .  .  (9) 

Both  the  probability (7) (x)  and  the  return  period  T  are  laid  off  on  auxilary 
horizontal  scales  on  the  probability  paper,  Eq.  7  relating  the  observed 
variate  x  to  the  reduced  largest  value  y  may  be  written  as 

X  =  -L.  y  +  U  (10) 

^n     n 

from  which  it  appears  that  there  is  a  linear  relationship  between  x  and  y. 

Thus,  theoretically,  the  plot  of  x  vs  y  should  be  a  straight  line  on  probability 
paper. 


(14)   Powell,  R.  IT. ,  "A  Simple  Method  of  Estimating  Flood  Frequencies," 
Civil  Eng.,  105-106,  1943. 


Benson,  '"'in  1950,  made  use  of  historical  data  in  flood-frequency 
analysis.  The  study  was  made  on  the  basis  of  annual  peak  discharges,  and 
data  were  plotted  on  probability  paper  based  on  Gurabel's  extreme  value 
theory. 

The  present  study  includes : 

1.  A  statistical  analysis  of  existing  peak  discharge  data  and  the 
determination  of  the  25-year  flood  for  32   gaged  small  watersheds 
by  Gumbel's  extreme  value  method „ 

2.  A  geomorphological  analysis  of  16  of  these  small  watersheds. 

3.  A  multiple  correlation  of  the  25-year  flood  and  the  geomorphological 
characteristics.  This  is  based  on  the  assumption  that  the  State  of 
Indiana  is  an  area  sufficiently  homogeneous,  so  that  variation  of 
storm  characteristics  and  of  physical  properties  of  soils  are  not 
significant  variables  compared  to  the  variables  describing  the 
geomorphological  characteristics.  The  watersheds  considered  in  this 
study  have  an  area  between  20  and  200  square  miles  approximately. 
These  watersheds  are  called  small  because  they  refer  to  the  smaller 
group  for  which  flow  records  are  reported  by  the  U.S.G.S.  They  are 
large  enough,  however,  so  that  the  land  use  or  the  type  of  vegetation 
is  not  an  important  variable,  and  may  be  disregarded.  This  method 
has  an  advantage  that  the  estimation  of  design  discharges  is  based  on 
the  discharge  records  themselves,  thus  eliminating  the  need  of  re- 
lating rainfall  and  runoff  by  means  of  a  number  of  variables  difficult 
to  evaluate. 


(15)  Benson,  M.  A.,  "Use  of  Historical  Data  in  Tlood-Frequency  Analysis," 
Amer.  Geophy,  Union  Trans,,  V.  21,   pp.  419-424,  June  1950* 
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Flood  frequency  analysis  for  small  water sheds  in  Indiana 

Thirty-two  watersheds  distributed  throughout  the  whole  state  were 
selected  to  study  for  flood  frequency  analysis.  Fig.  1  is  a  map  of 
Indiana  showing  the  thirty-two  watersheds,  and  Table  1  lists  the  names 
of  the  watersheds,  their  assigned  number  and  their  areas.  For  convenience, 
these  numbers  will  be  used  in  this  report  instead  of  the  name  of  the 
watersheds.  There  are  only  three  watersheds,  the  area  of  vhich  is  less 
than  50  square  miles,  12  watersheds  under  100  square  miles,  and  29  water- 
sheds under  200  snuare  miles  for  which  flow  records  are  available  from 
the  U.S.G.S.  All  of  these  were  included  in  this  freouency  analysis.   In 
order  to  have  a  good  distribution  of  the  watersheds  over  the  whole  state, 
three  additional  larger  watersheds  were  included.  A  bar  diagram  was 
plotted  in  Fig.  2  to  show  the  time  period  of  records  for  each  watershed. 

Instantaneous  annual  peak  flow  was  used  for  the  extreme  value  analysis 
instead  of  the  maximum  mean  daily  flow,  since  it  is  the  peak  discharge 
that  is  desired  for  hydraulic  structures  design.  Data  for  this  instantaneous 
annual  peak  discharge  were  obtained  from  the  U.S.G.S.  office  in  Indianapolis, 
Indiana.  The  data  are  analyzed  on  a  water  year  basis — that  is,  from 
October  to  September.  The  data  were  examined  for  man-made  changes,  such  as 
new  reservoirs,  soil  and  water  conservation  works,  and  changes  of  gaging 
sites.   It  was  found  that  no  correction  was  needed. 


Table  I  -  List  of  Watersheds,  Their  Area  and  Aligned  Number 


Watershed  Watershed  Area 

No.  Gaging  Station  (sq.  mi.) 


S-l  Bean  Blossom  Creek  at  Dolan,  Ind.  100 

3-2  Clifty  Creek  At  Hartsville,  Ind.  88.8 

S-3  North  Fork  Vernon  Fork  near  Butlerville,  Ind.     87.3 

S-4  Hart  Ditch  at  Munster,  Ind.  69.2 

S-5  Salt  Creek  near  McCool,  Ind.  78.7 

S-6  Little  Calumet  River  at  Porter,  Ind.  62.9 

3-7  Cedar  Creek  at  Auburn,  Ind.  93.0 

S-8  West  Creek  near  Schneider,  Ind.  46.3 

S-9  Iroquois  River  at  Rosebud,  Ind.  30.3 

S-10  Bice  Ditch  near  South  Marion,  Ind.  22.6 

S-ll  Big  Slough  Creek  near  Collegeville,  Ind.         84.1 

S-12  Carpenter  Creek  at  Egypt,  Ind.  48.1 

1  Tippecanoe  River  at  Oswego,  Ind.  115 

2  Mississinewa  River  near  Ridgeville,  Ind.  130 

3  Wildcat  Creek  at  Greentown,  Ind.  162 

4  Cicero  Creek  near  Arcadia,  Ind.  131 

5  Fall  Creek  near  Fortville,  Ind.  172 

6  Eagle  Creek  at  Indianapolis,  Ind.  179 

7  Young  Creek  near  Edinburg,  Ind.  109 

8  Blue  River  at  Carthage,  Ind.  187 

9  Sand  Creek  near  Brewersville,  Ind.  156 

10  North  Fork  Salt  Creek  near  Belmont,  Ind.  120 

11  Patoka  River  at  Jasper,  Ind.  257 

12  Busseron  Creek  near  Carlisle,  Ind.  228 

13  East  Fork  ,Thite  Water  River  at  Richmond,  Ind.    123 

14  Silver  Creek  near  Sellersburg,  Ind.  188 

15  Big  Indian  Creek  near  Corydon,  Ind.  129 

16  Kankakee  River  near  North  Liberty,  Ind.  152 

17  Singleton  Ditch  at  Schneider,  Ind.  122 

18  Deep  River  at  Lake  George  Outlet  at  Hobart,  Ind.  125 

19  Pigeon  Creek  at  Hogback  Lake  Outlet  near 

Angola,  Ind.  102 

20  Laughery  Creek  near  Farmers  Retreat,  Ind.  248 
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The  historical  floods  were  included  in  the  analysis  as  suggested  by 
Benson,' ■'-5)  Dalryraple , ^   '  in  order  to  extend  the  short  time  records 
available.  The  recurrence  interval  of  historical  floods  may  be  computed 
for  the  following  two  cases : 

1.  When  the  historical  flood  is  higher  than  those  for  the 
period  of  record.   In  this  case  the  recurrence  interval  of 
the  historical  flood  is  recorded  as  equal  to  one  plus  the 
number  of  years  of  the  period  for  which  it  is  the  greatest. 
For  example  the  1913  historical  flood  in  Patoka  River  at 
Jasper,  Indiana  (watershed  Fo.  11)  was  16,000  cfs.  The 
period  for  which  it  was  the  greatest  is  47  years  (1913  to 
1959).  According  to  Jq.   8-   the  corresponding  probability 
is  fj)  =  47/48,  and  the  recurrence  interval,  by  En,  9  is 

T  =  48  years. 

2.  When  the  historical  flood  is  known  to  be  the  highest  until  a 
greater  flood  occurs  during  the  period  of  record.   -An  example 
of  this  is  the  1913  historical  flood  in  Eagle  Creek  at 
Indianapolis,  Indiana  (watershed  Fo.  6).   It  was  19,000  cfs, 
and  was  the  second  highest  flood.  The  1957  flood  was  28,800 
cfs  which  was  larger  than  the  historical  one.  The  recurrence 
interval  of  these  floods  are  calculated  as  in  Table  2: 

Table  2   Calculation  of  Historical  Floods 


Floods  n      (£) 


28,800  cfs.  Max.  (1913—1959)  47    47/4*      48 

1^,000  cfs.  2nd.  Flood  (1913—1959)  47    46/48      24 


The  flood  frequency  analysis  based  on  Gumble's  extreme  value  theory 
was  done  for  all  32  selected  watersheds.  The  instantaneous  annual  peak 
discharges  were  plotted  against  the  reduced  value  y  on  the  probability 
paper,  and  a  straight  line  was  fitted  by  the  method  of  the  least  squares, 
except  for  watersheds  No,  4,  6,  14,  18,  and  20.  These  were  fitted  by 


(16)  Dalrymple,  Tate,  "Flood-Frequency  Analyses,"  U.S.G.S.  Water  Supply 
Paper,  1543-A,  I960. 


inspection  in  order  to  give  more  weight  to  the  lower  points  because  the 

upper  point  or  the  historical  flood  plotted  excessively  high  in  comparison 

with  the  straight  line  pattern  of  the  lower  points.  By  extending  these 

lines,  the  expected  floods  of  different  frequencies  were  obtained  and  listed 

in  Table  3 .  An  example  is  given  as  follows : 

Watershed  No.  S-6 

Gaging  Station:  Little  Calumet  River  at  Porter,  Indiana 

Computations  are  shown  in  Table  4,  and  plotted  results  are  as  shown  in 

Fig.  3. 

TABLE  3 

The  Predicted  Annual  Instantaneous  Peak  Discharge  from 
Flood  Frequency  Analvsis 


Watershed 

Fred 

icted  Annual 

Instantaneous 

Peak  Discharge 

No. 

25-yrs. 

50-yrs. 

75-yrs. 

100-yrs . 

cfs. 

cfs. 

cfs. 

cfs. 

S-l 

11800 

13800 

14700 

15900 

S-2 

12900 

15000 

16000 

17200 

S-3 

20500 

24000 

25400 

27000 

S-4 

3300 

3750 

3940 

4200 

S-5 

2950 

3420 

3600 

3860 

S-6 

1630 

1800 

4000 

4300 

S-7 

1630 

1800 

1880 

1980 

S-8 

2100 

2350 

2450 

2600 

S-9 

465 

520 

540 

570 

S-10 

835 

925 

960 

1010 

S-ll 

2440 

2770 

2900 

3080 

S-12 

4160 

4800 

5100 

5500 

1 

4160 

4800 

5100 

1110 

2 

140C0 

16500 

17400 

18800 

3 

7400 

8600 

9000 

9700 

4 

6800 

7800 

8300 

8900 

5 

7700 

8900 

9300 

10000 

6 

17000 

19700 

210C0 

22300 

7 

10500 

12400 

13000 

14000 

8 

10800 

12200 

12800 

13600 

9 

21500 

24000 

25200 

26800 

10 

19000 

22000 

23300 

25200 

11 

13500 

15600 

16500 

17700 

12 

8700 

9900 

10300 

11000 

13 

19800 

22800 

24200 

26000 

U 

13300 

15100 

15900 

17000 

15 

22300 

25700 

27300 

29400 

16 

1020 

1130 

1170 

1230 

17 

1270 

1370 

1410 

1470 

18 

4800 

5500 

5700 

6200 

19 

760 

860 

900 

960 

20 

28600 

33000 

35000 

3730© 
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TABLE  k 
Little  Calunet  River  at  Porter,  Indiana 


Year 

Instantaneous  Annual 
Peak  Discharge 
Q  cfs 

Rank 

o 
cfs 

<i> 

y 

1945 

2,440 

1 

490 

0.0625 

-1.01979 

1946 

715 

2 

521 

0.1250  -0.73210 

1947 

2,140 

3 

690 

0.1875 

-0. 51520 

1948 

1,960 

4 

715 

0.2500 

-0.32663 

1949 

690 

5 

848 

0.3125 

-0.15114 

1950 

1,720 

6 

1,060 

0.3750 

0.01936 

1951 

1,360 

" 

1,170 

0.4375 

0.19034 

1952 

1,060 

8 

1,360 

0.5000 

0.36651 

1953 

521 

9 

1,370 

0.5625 

0.55275 

1954 

1,170 

10 

1,420 

0.6250 

0.75501 

1955 

3,110 

11 

1,720 

0.6875 

0.98165 

1956 

1,370 

12 

1,960 

0.7500 

1.24590 

1957 

848 

13 

2,140 

0.8125 

1.57196 

1958 

490 

14 

2,440 

0.8750 

2.01342 

1959 

1,420 

15 

3,110 

0.9375 

2.74063 
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GEOIiORPHOLOGICAL  STUDY 

(17) 

Langbein  anH  others,   in -1947  irtdicatod  th-.t -river  TlnoTte  vere  'the 

results  of  many  causes.  One  of  the  primary  objectives  of  scientific 
hydrology  is  the  segregation  and  evaluation  of  the  causative  factors.  The 
climatic  factor  and  the  soil-vegetation  complex  are  variables  that  exercise 
their  principal  influence  on  the  volume  of  runoff.  The  topography  of 
drainage  basins  is  a  reasonpbly  permanent  characteristic  which  influences 
mainly  the  concentration  of  time  distribution  of  discharge  from  a  drainage 
basin.  The  topographic  characteristics  of  drainage  basin  cited  were:  area 
of  basin,  drainage  density,  area-distance  distribution,  length  of  basin, 
land  slope,  area-altitude  distribution,  and  area  of  water  surfaces. 

Strahler, (1°)  in  1952,  expressed  the  area-altitude  relation  by  a  hypso- 
metric analysis.  The  hypsometric  curve  relates  horizontal  cross-section 
area  of  a  drainage  basin  to  relative  elevation  above  basin  mouth.  By  the 
use  of  dimensionless  parameters,  the  curves  can  be  described  and  compared 
irrespective  of  true  scale.  The  area  under  the  curve  can  be  used  to  find 
total  land  mass  and  the  mean  relief  of  watersheds. 

Strahler, (19)  in  1957,  made  a  quantitative  analysis  of  watershed 
geomorphology  and  showed  that  the  linear  scale  measurements  include  length 
of  stream  channel  of  given  order,  drainage  density  and  relief.  Surface  and 
cross-sectional  area  of  basins  are  length  products.  Dimensionless  proper- 
ties include  stream  order  numbers,  stream  length  and  bifurcation  ratio, 


(17)  Langbein,  U.B.,  "Topographic  Characteristics  of  Drainage  Basin," 

U.S.G.S.  \rater  Supply  Paper,  968-c,  pp.  155,  1947. 

(18)  Strahler,  A.  N. ,  "Hypsometric  (Area- Altitude)  Analysis  of  Erosional 

Topography,"  Bui.  Geol.  Soc.  Amer.,  V.63,  pp.  1117-1142,  1952. 

(19)  Strahler,  A.  N. ,  "Quantitative  Analysis  of  Watershed  Geomorphology," 

Technical  Report  No.  13,  Dept.of  Geol.,  Columbia  University,  1957, 
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maximum  valley-side  slopes,  mean  slopes  of  watershed  surfaces,  channel 
gradients,  relief  ratio  and  hypsometric  curve  properties. 

Benson, W  in  1959,  cited  the  following  basin  characteristics:  drainage 
area,  channel  slope,  land  slope,  tributary  channel  slope,  watershed  shape 
factor,  mean  elevation,  percentage  of  lakes,  swamps  and  reservoirs,  and 
drainage  density.  The  shape  factor  of  the  watershed  was  represented  by: 

A   ,  equivalent  of  basin  length  divided  by  basin  width. 


L 

A   ,  basin  length  divided  by  drainage  ar°a, 

'S  aL  ,  the  summation  of  small  subdivisions  of  the  drainage  area, 

each  multiplied  by  the  distance  of  travel  to  the  gaging  point. 
In  this  study,  the  first  geomorphological  factor  considered  is  the 
drainage  area.   It  is  the  projected  area  of  watershed,  also  called  the 
catchment  area.  Obviously,  the  bigger  the  area  of  catchment  the  larger  the 
amount  of  runoff.  However,  the  rate  of  runoff  is  largely  dependent  on  the 
slope  of  the  land,  the  drainage  density,  and  the  slope  of  the  streams.  A 
steep  land  slope  with  high  drainage  density  and  larger  slope  of  stream  will 
give  a  higher  rate  of  runoff  than  those  which  have  smaller  slope  and  lower 
drainage  density.  The  shape  of  the  watershed  is  also  a  factor  affecting  the 
runoff,  becuase  it  affects  the  time  of  concentration.  With  a  three-dimen- 
sional concept  in  mind,  the  principal  geomorphological  factors  which  can  af- 
fect the  amount  of  peak  discharge  are  listed  as  following. 

1.  The  drainage  area 

2.  The  drainage  density 

3.  The  land  slope 

4.  The  main  stream  slope 

5.  The  watershed  slope 

13 


Evaluation  of  geomorphological  factors 

1.  Drainage  Area 

Area  of  the  watershed  is  directly  measured  from  topographic  maps 
with  a  planimeter.   It  is  expressed  in  square  miles. 

2.  Drainage  Density 

Drainage  density  is  defined  as  the  total  length  of  streams  in  the 
watershed  divided  by  its  total  area,  that  is  the  length  of  streams  per 
unit  area  of  the  watershed.   It  can  be  expressed  as  follows: 

Drainage  Density  (See  Figure  U) 


D 


£L 


A  ■  (11) 

IJhere^L  is  the  total  length  of  streams,  and  A  is  drainage  area. 
The  length  of  stream  is  expressed  in  miles  and  is  obtained  from  drain- 
age maps.  Since  the  area  is  expressed  in  souare  miles,  the  drainage 
density  can  be  expressed  as  miles  per  souare  miles. 
3.  Land  Slope 

Since  the  land  slope  is  changing  from  place  to  piece  in  a  water- 
shed, it  is  difficult  to  find  a  ouantitative  value  to  represent  the 
land  slope  of  a  whole  watershed.  A  new  parameter  introduced  here  to 
replace  the  land  slope  is  the  mean  relief  of  land.  The  mean  relief  is 
defined  as  the  total  volume  of  land  mass  above  the  outlet  of  a  water- 
shed divided  by  its  projected  area.  This  can  be  evaluated  auantita- 
tively  by  using  the  so-called  hypsometric  curve  aeveloped  by  Langbein 
and  by  others  which  give  a  dimensionless  relationship  between  the 
horizontal  cross-sectional  drainage  basin  area  and  the  elevation. 
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Fig.  5  shows  a  watershed  and  its  horizontal  projection.  The  maximum 
elevation  H'  may  be  obtained  from  topographic  maps,  and  the  cross- 
sectional  area  A  of  the  watershed  can  be  measured  by  planimeter. 
Similarly,  the  projected  area  above  any  height  h  may  be  obtained  from 
topographic  maps,.  The  dimensionless  plot  of  the  relative  area  a/A  against 
the  relative  height  h/H'  is  called  a  hypsometric  curve,  and  its  general 
aspect  is  as  shown  in  Fig.  6e   From  the  hypsometric  curve,  the  total  volume 
of  land  mass  and  the  mean  relief  can  be  calculated.  Since  the  area  under 
the  hypsometric  curve  can  be  easily  measured,  this  can  be  expressed  as, 


which  is  rewritten  as 


J1!  *<|,  >-«!  (12) 


1  r1 

IhTI  adh   =    c*1  (13) 


Hence  the  total  volume  of  land  mass   is 

V  =    f1     adh     =  AH'oc1  (U) 

Jo 

and  the  mean  relief  is 

V  i 

H=^-cvH'  (15) 

It  is  thus  seen  that  the  mean  relief  is  equal  to  the  product  of  the 

area  under  the  hypsometric  curve  and  the  maximum  height  over  the  outlet. 

4.  Main  stream  slope 

The  slope  of  the  main  stream  can  be  obtained  from  the  topographic 

map  at  several  points  along  the  length  of  the  stream.  Usually,  the 

upper  stream  reaches  are  steeper,  and  the  downstream  reaches  are  flatter. 

The  mean  slope  is  calculated  by  means  of  the  formula  introduced  by 


15 


(20) 

Tavlor  and  Schwa rz: 


S  = 


(16) 


,1*1+1    +   1 
(sl/2   gl/2    gl/2   •••gl/2) 


°n 

where  n  is  numbers  of  equal  reaches,  S.  ,  Sp,  S  , and  S  ,  are  the  slopes 

2/3 

of  each  reach.   It  is  based  on  the  assumption  that  the  quantity    R 

n 
which  appears  in  the  Manning  formula,  is  the  same  in  all  reaches. 

5.  Watershed  shape  factor 

The  watershed,  shape  factor  is  the  ratio  of  Lhe  main  stream  length  to 

the  diameter  of  a  imaginarv  circular  watershed,  of  equal  area.  Figure  7 

shows  two  watersheds  with  the  same  area  but  different  shape,  the  actual 

irregularly/  shaped  watershed  and  the  imaginarv  circular  one.  A0  is  measured 

along  the  main  stream  up  to  the  watershed  boundary  line,  and  A'O1  is  calculated 

from  the  known  area  of  the  actual  watershed 


a'C  =2    /  Area  of  actual  watershed (17) 

pi  E 

The  shape  factor  is  determined  as 

f  =  -£0.  (18) 

A'O' 

GeomorpholoFical  factors  of  small  watersheds  in  Ind iana 

The  followinr  Table  5  shows  quantitatively  the  geomorphological  factors 

of  16  small  watersheds  in  Indiana,  the  only  ones  for  which  topographic  maps 

were  available.   The  stream  length  and  stream  density  were  obtained  from 

(21) 
the  drainage  maps  of  Indiana. 


20 

Taylor,    <v.   A.   and  Schwarz,   H.   E.,    "Unit-ffydro graph  Lag  and  Peak  Flow 

Related   to  Basin  Characteristics,"  Trir«.    A.G.U.,   V.   33,   pp.    235-246, 
1952. 

21 

Purdue  University,  "A.tlas  of  County  Drainage  MapS   Indiana,"  Joint 

Highway  Research  Project,  Engineering  Bulletin  Extension  Series. 
No.  97',  July  1959- 
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t  .kle  5 

atershed   Characteristics  and  25-Year  Annual  Instantaneous 
peak  Runoff  of  16  Watersheds  in  Tndi'-r.a 


'  ytershed 

No. 


25-Year  innual 
Instantaneous 
Peak  Runoff 
Q-cfs. 


atershed 


Characteristics 


vrea 

.  so.  mi 


Mean 
Relief 
H.  ft. 


Drainage 

Density 
d  mi/sqlri 


Shape 
Factor 
f 


Main 
Strean 

Slope 
3x  10' 


-4 


-u 


S-l 

11,-00 

100 

216 

10.66 

2.63 

9.84  xlO 

3-2 

12,900 

88.8 

270 

7.3C 

3.00 

20. 38 

S-5 

2,950 

78.7 

101 

6.57 

1.75 

9.05 

S-6 

3,300 

62.9 

110 

8.00 

1.12 

21.10 

3-7 

1,630 

93.0 

79 

5.10 

1.47 

8.29 

1 

cvq 

115 

65.4 

3.35 

1.41 

2.64 

6 

17,000 

179 

195.2 

7.88 

2.2 

13.40 

7 

10,500 

109 

86 

7.02 

1.94 

10.39 

9 

21 , SCO 

156 

250 

9.°6 

2.85 

10.68 

10 

19,000 

120 

237 

11.20 

2.18 

9*90 

11 

13,500 

2r>7 

181.5 

13.95 

2.87 

2;95 

12 

:1,7C0 

228 

99.8 

13.20 

1.93 

5.43 

14 

13,300 

188 

195.8 

10.47 

1.35 

6.21 

15 

22,300 

129 

231 

8.70 

2.56 

10.16 

18 

4,800 

125 

•84-.  7 

4.50 

1.91 

6U35 

L9 

760 

102 

66.1 

3.16 

1.93 

7.93 
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Multiple  correlation  of  peak  discharge  and  geomorphological  characteristics 

A  multiple  correlation  was  derived  between  the  25-year  flood  and  the 
geomorphological  variables  considered  above.  Such  a  correlation  presumes 
that  the  area  of  application  is  meterologically  and  geologically  homo- 
geneous, otherwise  the  state  should  be  divided  into  different  zones  in 
which  homogeneous  conditions  exist. 

Since  the  State  of  Indiana  is  relatively  flat,  the  orographic 
precipitation  is  not  a  factor  in  the  larger  storms.  Convectional  pre- 
cipitation resulting  from  most  thunderstorms  usually  has  a  duration  which 
is  less  than  the  time  of  concentration  for  watersheds  of  more  than  20 
souare  miles.  Thus  frontal  and  cyclonic  precipitation  are  the  cause  of 
the  large  storms  producing  peak  runoff.  Therefore,  it  could  be  assumed 
that  the  climatolo^ic';!  condition  is  homogeneous  over  the  State  of  Indiana. 
That  is,  the  probability  of  being  subjected  to  a  storm  of  a  given  frequency 
is  almost  enual  for  all  watersheds  in  Indiana. 

By  studying  the  soil  map'^2)  0f  the  State  of  Indiana,  it  was  found  that 
the  permeability  varies  from  soil  to  soil.   Although  this  would,  in  theory, 
disprove  the  assumption  that  the  geological  conditions  are  homogeneous 
throughout  the  state,  one  additional  multiple  correlation  including  the 
maximum  intake  rate^  -'    of  the  soil  as  a  variable  indicated  that  it  was  not 
significant  compared  to  the  geomorphological  variables. 


(22)  Belcher,  D.J.,  Gregg,  L.E.,  and  Woods,  K.B.,  "The  Formation,  Dis- 

tribution and  Engineering  Characteristics  of  Soils,"  J »int  High- 
way Research  Project,  The  State  Highway  Commission  of  Indiana, 
and  Purdue  University.  January,  1943. 

(23)  Purdue  Universitj^,  "\  guide  for  Designing  Sprinkler  Irrigation  Systems 

in  Indiana,"  Agricultural  Extension  Service  in  Cooperation  with 
Soil  Conservation  Service,  U.S.u.jl,  1955. 
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(24) 

Multiple  correlation    is  a  statistical  method  to  find  the  relation- 
ship between  one  dependent  variable  and  a  number  of  independent  variables. 
If  a  linear  relationship  exists,  the  method  of  fitting  is  to  make  the 
sum  of  the  squares  of  the  deviations  of  actual  observations  from  the 
theoretical  linear  relation  a  minimum.  This  is  called  the  method  of 
least  sciuares. 

The  freneral  formula  for  multiple  correlation  is 

y  =  a  +  b1x1  +  b^  +  b^  (19) 

where  y  is   the  dependent  variable,   Xp   x2,   x-,, x,    are  the   independent 

variables,   and  a,   b-,,   b2,   b   , b   ,    are   constants  obtained  by    solvinr  the 

following    si  ultroneous  equations: 

bx  S(xx)2  +  b2  S(Xlx2)   +  b3  S(x1x3)    +  bk  SCr^)   =  S(xiy) 

b±  S(x2Xl)   +  b2  S(x2)    +  b3  5   (x2x3) +  \  3(x2xk)   "  S   (x2y) 

h,_  s(x3xi)  +  b2  s^)  +  b3  s(x2)  .......; ;.+b^  s(x^}  m  s{^7) 

b1  SCx^)   +  b2  5(xkx2)   +  b3  S{x^) +bR  3(xR)   =  S(xky) 

and 

a  =  y  -  b^     -  b  x       -  b^    -b^k  (20) 

where  S(xy)   is  a  svmbol  which  means  J   xy  -  "—  x^»- v 

n 

so  that 

v2 


n 
S   (xxx2)   =  2Ix1   x2  -     £X1^X2 

5   (Xlxk)   =-x1xk  -  _ 


2k 

Bennett,   C.   A.   and  Franklin,   N.  L.,    "Statistical  Analvsis  in  Chemistry 
and  the  Chemical  Industry"  John  ,rilev  &  Sons,    Inc.   1954. 
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The  general  formula  for  estimate  of  variance  is 

S*   >-.vl  ■ =  s^y2)  -^   S^-^   "b2  S(x,y) -\  S(xky) 

?,   xl>  x2>  x3  xk  „   ,. 

n  -  (k  +  1) 

(21) 
where  k  is  number  of  independent  variables. 

This  serves  as  a  measure  of  the  degree  of  correlation.  Generally 

the  smaller  the  variance,  the  better  the  correlation. 

The  following  Ta^le  6  shows  the  results  of  correlation  between  the 

25-""ear  instantaneous  peak  discharge  and  the  geomorphological  factors 

of  small  watersheds  in  Indiana.   The  regression  formulas  in  exponential 

type  as  shown  in  Table  6  were  obtained  by  means  of  a  logarithmic 

transformation  and  the  formulas  (19)  and  (20).  The  standard  deviation 

were  calculated  by  means  of  Eq.  (21). 

Table  6  -  hultiple  Correlation  between  Predicted  25-years 
Instantaneous  Peak  i)isch<i.rge  and  G-o^orphological  factors 

Discharge  ^Gecncrph<.»lug±ca_     •Kogregsion  for;  .ul-             '    standard*-    ~ 
—. 1  ,-ctors D^viat.  inn 

0.847!  1.9?45 
Q       A,  f      ~  =  30.29  A      f  0.394 

.  .  792  1.5*29 
A,  3       0  =  0.0005716  i      S  0.305 

0.1994   1.9847 
a,  D      q  -  47.13  A      D  0.296 

0.7220  1.7396 
A,  H      Q  =  0.03993  A      H  0.234 

0.7066  0.1701  O.IICG 
0        A,  H,  t    0  =  0.48  A      H      f  0.243 

0.3036  1.2732  0."601 
L,  H,  D    r    =  0.31-32   S      H      D  0.213 

1.4623  1.3035  0.6938 
k,   B,  S    0  =  0.0022  S      H      S  0.211 

1.3915  1.0825  0.7114  0.6339 
0        A,  H,  S,  f.Q  =  C.0057  A     H      S      f  0.217 

1.0513  0.9672  0.5<390  0.7417 
A,  H,  3,  D  0  =  0.02337  \  H      3      D  0.199 

0.9715  0.7344  0.5901 
Q        A,  P.,  S,  D,  f  Q  =  0.05363  A      H      3 

0.3234  C.4160  0.190 

D      f 

*  in  Log-unit. 
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A  whole  set  of  multiple  correlations was  listed  in  Table  6.  By  studying 

the  standard  deviations,  it  is  e-.sv  to  find  the  significance  level  of 

reomorpholoricc-.l  factors  v.'hicb  affect  the  discharge.   If  A  is  assumed 

to  be  the  first  in  significance  to  influence  the  discharge,  then  the 

H  will  be  the  second  one,  and  S,  D,  f  will  be  the  third,  fourth  and  last 

respectively-.  The  last  regression  formula  in  Table  6} 
0.9715  0.781A.   0.82U  O.Z;l6C  0.5901- 
q  =  0.05363  4      F      D      f      3  (22) 

with  the  least  standard  deviation  is  of  course  the  "'-est  expression  obtained 

from  multiple  correlation.   If  all  the  geornorphological  factors  with  their 

powers  together  are  combined  as  a  "basin  characteristic",  B, 
C.9715  0.  e,irL   0.;VO34  O.'loO  0.5901 
B  =  X  F      D      f      S  (23) 

then,  the  regression  formula  can  be  expressed  as  H   =  cB,  where  c  is  a 

constant.  This  is  the  equation  of  a  straight  line  on  log-log  paper 

as  shown  in  Figure  8.  A  95"  confidence  interval  was  calculated  and  plotted. 

The  probability  of  anv  obsei /ation  falling  between  these  intervals  is  95c"*. 

To  test  the  assumption  of  geological  horn ronoity  an  additional  multiple 

co -relation  was  mace  including  a  variable  representing  soil  characteristics. 

(23) 
The  maximum  intake  rate    of  the  soil  expressed  in  inches  per  hour  was 

selected  as  the  new  variable  in  addition  to  the  five  geornorphological 
variables  used  before.  The  six  variable  multiple  correlation  which 
would  take  care  of  the  nonbonogeneons  condition  of  the  watershed  geology 
did  not  show  much  improvement  in  the  degree  of  correlation  to  the  peak  dis- 
charge, having  a  standard  deviation  somewhat  larger  than  that  for  Eq.  (22). 
This  is  because  the  variation  in  the  intake  rate  is  small  for  the  watersheds 
studied.  Until  additional  data  on  runoff  from  small  watersheds  becomes 
available,  En.  (22)  is  the  best  relationship  that  can  be  derived  between 
peak  flow  and  physiographic  factors. 
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Construction  of  correlation  cnart  for  predicting  future  flood 

The  regression  formula,  Eq.  (22)  or  its  graphical  representation  in 
Figure  2  may  be  used  to  calculate  the  flood  discbarge  if  the  geomorphological 
factors  are  known.   \  nomographic  representation  of  the  regression  formula, 
Eq.  (22), was  prepared.   From  this  it  is  possible  to  obtain  directly  the 
flood  discharge  from  the  knoxm  geomorphological  factors. 

The  regression  formula  has  the  form 


Letting 


then 


0.  =  KAa  Bb  Gc  Dd  Ee 


M  =  Aa  Bb 


Rj  =  Aa  Bb  GC  =  M  CC 


M  =  Aa  Pb  C°  Dd  =  M9  Dd 
3  2 


0  =K/aBb  CC  Dd  E6  =KlI3Ee 


(24) 


The  above  eouations  may  be  plotted  with  3,  C,  D,  E  as  parameters 
and  b,  c,  d,  e  are  knovn exponents  of  the  regression  formula  as  shown 
in  Figure  9.   The  diagrams  of  Fig.  9  may  be  combined;  the  "A  "  scale 
is  changed  from  "A"  scale  so  that  the  area  a   can  be  introduced  directly 
into  the  chart  along  with  the  other  factors  necessary  to  find  the  flood 
discharge,  as  shown  in  Figure  10. 

The  above  procedure  was  used  to  prepare  the  nomograms  of  Figures  11 
and  12.  These  may  be  used  instead  of  Eq.  22  or  of  Figure  %   to  estimate 
the  25-year  peak  discharge  prom  the  geomorphological  characteristics.   Fig.  11 
was  prepared  for  •   watersheds  less  than  100  square  miles  and  in  Fig.  12, 
the  correlation  is  extrapolated  to  cover  watersheds  up  to  300  square 
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miles.  Two  examples  are  presented  to  show  the  use  of  the  charts: 
(a)  Watershed  To. 
Watershed  area 
Mean  relief 
Drainage  density 
Sh^pe  factor 
Main  Stream  slope 
Following  the  dotted  line  in  Pip.  11,  the  flood  discharge  is  read 
from  the  chart  direct!^  as  3,'750  cfs.,  vhile  the  flood  predicted 
by  the  frequency  study  is  3,30C  cfs. 


5-6 

(0 

62.9  so.  mi. 

(H) 

110  ft. 

(D) 

8  mi./sq.  mi. 

(f) 

-4 
(3)  x  10 

1  12 

• 

21.10  x  10 

(b) 

Watershed  To. 
Watershed  area 
Mean  relief 

(-0 
C) 

3-1 

100  sq.  mi. 

216  ft. 

Drains  re  ^ensit-"- 

(D) 

10.66  mi./sq. 

mi. 

Shape  factor 

(f) 

-u 

2.63 

Main  stream  slope 

(3) 

x  10 

9.84  x  10 

Following  the  dotted  line  ir.  rig.  12,  the  "lood  discharge  reading 
is  12,' 50  crs.,  while  the  ."lood  predicted  by  the  frequency  study 
is  11,800  cfs. 
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L  simple  approximate  estimation  of  peak  discharge  'or  small  watersheds 

in  Indiana. 

As  shown  in  the  previous  paragraphs,  equation  (22)  is  the  best  one 

among  the  10  formulas  listed  in  Table  6  "or  predicting  the  peak  discharge 

for  sm:  11  watersheds  in  Indiana  since  it  has  the  least  standard,  deviation. 

Although  a  graphical  representation  of  the  correlation  formula,  Eq.  22 

is  given,  its  use  is  time  consuming  because  of  the  tedious  work  required 

for  the  deteiTninati^n  of  the  :"ive  geomo  phological  factors  which  enter 

into  the  formula  or  in  the  correlation  charts  as  independent  variables. 

In  particular  the  determination  of  the  drainage  density  D  and  of  the  mean 

relief  F  are  time  consuming,   ""or  the  multiple  correlations  listed  in 

Table  6,  it  appears  that  the  formula  containing  i,  F  and  S  has  a  standard 

deviation  of  0.211  which  is  close  to  the  standard  deviation  of  0.190  for 

Eq.  22  and  mav  therefore  be  use^  as  a  good  approximation  in  practical 

design.  This  formula  eliminates  the  neeo1  of  calculating  the  drainage 

densitv  D,  and  the  ivaters^ed  shape  factor  f  which  is  the  least  significant 

variable  is  omitted .  Tbe  formula 

1.1,621       1.3035    0.6938 
«*  =  0.0022  \  F        3  (25) 

still  requires  the  determination  of  the  mean  relief  F  of  the  watershed. 

however,  there  is  an  approximate  way  of  determining  F  by  estimating  ^ , 

the  area  under  the  hypsometric  curve.   The  average  hypsometric  curves 

have  been  plotted  "or  small  watersheds  in  Indiana  in  Fig.  13,  from  which 

it  appears  that  the  o£  -  values  vary  from  0.4  to  0.3.   From  the  pattern  of 

the  hypsometric  curves  it  is  possible  to  estimate  o<  by  calculating  only 

one  value  of  a/'v  for  a  certain  elevation  ratio,  say  h/F'  =  0.5.  Thus 

two  area  measurements  are  needed:   the  total  area  of  the  watershed  v,  and 

the  area  at  an  elevation  half  wav  between  the  maximum  elevation  and  the 

mouth.  The  following  Table  7  gi  res  the  CK.   values  "or  corresponding  a/n: 
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Table   7     -     The  O^  values   corresponding  to  a/A 
for  h/F' r  0.5 


iA 


ex 


C.3  0.40 

0.4  0.45 

0.5  0.50 

0.6  0.55 

0.7  0.60 

0.8  0.'70 

0.9  0.c>0 


For  a  rourh  estimate  of  q^  ,    it  ria-r  be  assumed  that  it  varies  between  0.4 
to  0.  '  for  F-shape  valleys  end  'rom  '  .6  to  0.  7  for  '/-shape  valleys  in  Indiana,   ith 
the  vnlue  of   determined  from  Table  7,  the  nesn  relief  F  c::n  be  easily  calculated 
from  !  q^  En.  15,  and  the  peak  discharge  c:.n  be  obtained  by  Eq.  25  . 

Thus,  the  peal  discharge  cm  be  clcul-ted    three  methods.  The  first 
is  the  simplest;  c<  is  obtained  from  Table  7,  from  which  by  Eq.  15  and  Eq.  25 
a  first  approximation  of  the  peak  discharge  is  obtained.   In  the  second  method 
F  is  determined  by  plotting  the  hypsometric  curve  and  Eq.  25  is  used  as  a 
second  approximation  of  the  peal-  discharge.  The  third  method  requires  the 
evaluation  of  the  five  p;eomorphological  factors,  after  which  the  peak  discharge 
is  calculated  by  En.  22  or  by  means  of  figures  11  and  12. 

An  example  is  riven  l^elow  and  solved  fr"  the  three  methods. 

''atershed  No.   3-1 

First  Method: 

geornorphological  factors 

A  =  100  square  miles 
S  =  9.^4 
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H»  =400  feet 

at  h/H'  =  0.5,  a/A  =  0.573,  from  Table  7,  cl=  °.53 
H  =  400  x  0.53  =  212  feet 
-':;•-  Jq-,.  25,  the  25-year  peak  discharge  is 
Q  =  9,280  cfs. 
Second  Method: 

geomorphological  factors 

A  =  100  souare  miles 
S  =  9.84 
K  =  216  feet 
by  .V'.  25,  the  25-year  peak  discharge  is 
Q  =  9,^80  cfs. 
Third  Method: 

The  25-year  peak  discharge  as  found  by  using  the  regression  formula, 
Eq.  22  or     Figures  11  and  12  as  shown  in  the  Previous  cdar.   le  relating  to  those 
charts,  is  n   =  12,450  cfs. 

while  the  25-year  peak  discharge  estimated  from  the  frequency  study 
was  11,800  cfs. 

Comparing  the  results  obtained  by  the  three  methods  to  the  peak  dis- 
charge obtained  from  the  frenuencv  study,  it  appears  that  the  percentages  of 
error  are  respectively  -21.4$,  -19.7',  and  +5.5^.    It  should  be  remembered 
however  that  these  percentage  will  vary  from  one  watershed  to  the  next  as 
methods  1  and  2  are  calculated  by  means  of  one  regression  formula  whereas 
method  3  is  obtained  by  a  different  regression  formula. 


26 


Relationship  betv/een  25-vear  peak  discharge  and  the  peak  discharge  for 
other  freouencies. 

In  the  preceding  paragraphs,  the  peak  discharges  from  small  watersheds 
were  obtained  for  a  recurrence  interval  of  25-years  which  was  based  on  the 
average  life  of  small  highway  drainage  structures.  However,  it  may  be 
desirable  to  estimate  the  peak  discharge  for  other  return  periods  so  that 
the  design  engineer  may  have  a  greater  freedom  of  choice.  Hence  the  rela- 
tionship between  the  peak  discharge  for  any  freouency  and  the  25-year  peak 
discharge  was  derived  for  small  watersheds  in  Indiana.  The  thooretic-il  re- 
lationship is  based  on  Gumbel's  extreme  value  theory.  Figure  14  shows  two 
theoretical  straight  lines  for  any  two  watersheds.  The  differences  be- 
tween the  25-year  peak  discharge  and  the  n-year  peak  discharge   for 
the  two  watersheds  are  obtained  from  En.  (10),  and  are: 


1 


.&y  (26) 


C*2 

where  l/^,  1/fy,,  are  the  slopes  of  the  straight  lines,  and  the  increment  a  yof 
the  reduced  variate  corresponds   to  the  selected  increment  of  frequency. 
For  a  fixed  freouency  n  the  increment  ^y  is  a  constant.  A  general  form 
thus  can  be  written  for  all  the  watersheds  as 

(27) 


or 


=  c<o   (1  +  -***-  ) 

.-••.  )-25 

thus      n  =  o   (i  +   *Y     )  (28) 

'  ^25 
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an 

—  . 

i 

&y 

*    n 

-« 

n25 

B&y 

*    n 

=y. 

n25 

+  *y 

Taking  logarithms  on  both  sides 

log  Q     =   log  C        +  log   (1   +        A7      )  (?a) 

n  <^5  C<c25 

Eq.   29  is  a  linear  on  log- log  paper,   if  the  last  term  is  constant. 

An  examination  of  the  last  term  shows  firstly  that  the       Ay  is  sman 

o<025 
compared  to   "1"   and  secondly  that  the  value  ofOfQ2r  for  most  of  the   small 

watersheds  in  Indiana  varies  from  4  to  6  and  an  average  value  of  5   can  be 

used  forOfQ       for  small  watersheds  in  Indiana.     Then,  Equation  29  becomes 

l*g  Qn  =  log  Q25   +  log   (1  +  |y_)  (30) 

Values- of  Ay  for  several  aesign  frequencies  are  given  in  Table  8.  Values 
of  Ay  for  other  frequencies  may  b^  obtained  frotB  Fig.  3. 
Table  8.  Values  of  Ay  for  different  frequencies. 


n 

±y 

IC 

-o.93 

bO 

0.70 

75 

1.00 

100 

1.40 

Fig.  15  is  a  plot  of  Eq.  30  on  logarithmic  paper,  giving  the 
relationship  between  the  25-year  peak  discharge  and  the  peak  flow  for  frequencies 
of  10,  50,  75  and  100  -  veers,  oirfh  po.-Jc   discharges  may  be  determined 
from  Fig.  15  or  from  Eq.  (30)  if  the  25-year  discharge  has  previously  been 
determined  from  rFifuresll  and  12  or  from  Eq.  (22). 
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DI3CU33I0N  MD  CONCLUSIONS 

1.  The  paper  considers  watersheds  between  20  and  200  square  miles  in  area 
in  the  state  of  Indiana,  .-ill  available  past  observations  of  annual  peak 
discharge  were  plotted  on  a  probability  paper  using  Gumbel's  extreme 
value  theory,  as  there  is  a  linear  relation  between  the  observations 
and  the  reduced  largest  value  y,  the  best  fit  straight  line  was  then 
obtained.  Future  floods  with  different  frequencies  were  obtained  by 
extending  the  straight  line.   Table  3  gives  the  predicted  flood  of  25, 
50 j  75  and  100  years  frequency  for  3^  gaged  ./atersheds  in  Indiana. 

2.  The  quantitative  study  of  geomorphological  factors  of  small  water- 
sheds in  Indiana  combined  with  the  use  of  multiple  correlation  techniques 
gives  an  indirect  determination  of  peak  discharge.   This  is  based  on 

the  assumption  that  the  climatological  and  geological  conditions  are 
homogeneous  throughout  the  state.  Hence,  the  geomorphological 
characteristics  are  the  dominant  factors  which  affect  the  peak  dis- 
charge from  small  watersheds.   The  geomorphological  factors  considered 
significant  are:  the  watershed  area,  the  drainage  density,  the  mean 
relief  of  watershed,  the  main  stream  slope,  and  the  shape  factor  of 
the  watershed. 

3-   Correlation  charts  (Figs,  11  and  12)  were  prepared  to  obtain  the  25-year 
peak  discharge  directly  from  the  five  watershed  characteristics,  for 
areas  up  to  300  square  miles,   as  shown  in  the  previous  examples,  the 
design  engineers  may  use  these  design  charts  to  estimate  very  rapidly 
the  25-year  peak  discharge  .lith  good  accuracy.   The  peak  discharge  for 
other  frequencies  may  be  obtained  from  Fig.  15* 

4.  A  simple  formula  which  contains  only  three  geomorphological  factors , 

h,   H,  and  S  was  introduced  as  a  first-approximation.  H  is  determined 

by  the  ^  -  value  which  is  estimated  from  the  average  hypsometric  curves 
for  small  watersheds  in  Indiana. 
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Appendix  -  Notations 

A  Area  of  watershed,  sq.  miles 

a  Area  of  watershed  for  a  given  elevation  h,  sq.  miles 

D  Drainage  density  (see.  eq-.  11),  milcs/sq.  miles 

f  Watershed  shape  factor  (see  eq«  18) 

H*  Maximum  height  of  watershed  above  the  mouth,  ft. 

H  Mean  relief  of  watershed  (see  eq.  15),  ft. 

h  The  height  above  the  mouth  of  any  given  point  on  the  watershed,  ft, 

i_L  Total  length  of  streams  on  the  watershed,  miles 

m  Rank  of  the  yearly  maximum  in  increasing  order 

u  Uumbera  -of  independent  observations 

Q  Peak  discharge,  cubic  feet  per  second 

Q25  25-year  peak  discharpe,  cubic  feet  per  second 

°n  n-year  peak  discharge,  cubic  feet  per  second 

S  Mean  slope  of  main  stream  (see  ea  .  16) 

T  The  return  period,  ye^rs 

u  A  certain  expected  largest  value  having  the  return  period  n 
(see  eq  .  7) 

V  Total  volume  of  land  mass  of  watershed,  (see  eq1 .  14) 

x  Any  independent  observation  of  the  peak  discharge 

y  Reduced  largest  value  (see  eq~,  7) 

o(x)  Probability  of  x  being  the  largest  among  n  independent  ob- 
servations (see  eq  .   6  and  8) 

0>  Area  under  hypsometric  curve  (see  eq...  12) 

kJ  An  extreme  value  parameter  (see  en.  7) 
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Fig  2      PERIOD    OF  RECORD    OF    INSTANTANEOUS 
ANNUAL    PEAKS     AT    GAGING    STATIONS. 

■- 
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PERIOD   OF   RECORD 

* 

n 

1900 
1910 

1920 

1930 

1940 

1945 

1950 

1955 
I960 

S-  1 

14 

S-2 

12 

S-3 

18 

S-4 

17 

S-5 

15 

S-6 

15 

S-7 

17 

S-8 

9 

S-9 

II 

S-IO 

II 

S-M 

10 

S-12 

10 

1 

_ 

II 

2 

13 

3 

16 

4 

5 

5 

18 

6 

1 

21 

7 

17 

e 

M 

10 

9 

12 

10 

| 

14 

ii 

1 

1 

12 

12 

16 

13 

I 

10 

14 

5 

15 

17 

16 

9 

17 

II 

18 

13 

19 

14 

20 

19 

*    n=    Length    of   Record    in    Years 
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FIG.   3 
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